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Abstract 

We study a one-dimcnsional discrete analog of the von Karman 
flow, widely investigated in turbulence. A lattice of anharmonic oscil- 
lators is excited by both ends in order to create a large scale structure 
in a highly nonlinear medium, in the presence of a dissipative term 
similar to the viscous term in a fluid. This system shows a striking 
similarity with a turbulent flow both at local and global scales. The 
properties of the nonlinear excitations of the lattice provide a partial 
understanding of this behavior. 

1 Introduction 

Although studies of turbulence have been carried out for decades, they go on 
now-days because the problem is very hard and challenging but new ideas 
have emerged in the last few years from the investigations of statistical prop- 
erties of turbulence, both at a local or a global scale. First, studies of local 
longitudinal velocity differences measured at distance r showed that their 
probability densities are very well described by a theoretical model based on 
nonextensive statistical mechanics Q, which, in the context of turbulence, 
can be justified by a phenomenological model based on dynamical founda- 
tions Second, interesting properties showed up in the analysis of the 
statistical properties of global quantities such as the power consumption of 
the confined turbulent Von Karman fiow Q and they culminated with the 
observation of the universality of rare fluctuations in turbulence and critical 
phenomena Q. 

These recent developments are based on an analysis of experimental 
results [||, 0]. Indeed only such data contain the full features of turbulence. 
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but for an understanding of the results it would be desirable to have a simple 
model, showing the basic properties of turbulence, that would be amenable 
to very detailed studies. The search for simple models of turbulence lead to 
the introduction of the so-called shell models Q which describe the energy 
exchange mechanism between various scales. However shell models are not 
suitable to study properties that depend on space or to analyze how the 
dynamics can lead to the observed statistics. 

This is why it might be interesting to introduce models of turbulence 
that are sufficiently simple to allow a detailed study, and nevertheless pre- 
serve a spatial structure. In this letter we introduce such a model derived 
from studies on one-dimensional nonlinear lattices. These systems show 
various features which are also found in turbulence, and, in particular an 
exchange of energy between various spatial scales which tends to cause a flow 
of energy from large scales (long wavelengths) to smaller scales, through the 
mechanism of modulational instability . The idea to use one-dimensional 
systems as simple models for turbulence is not new and continuous one- 
dimensional media have been introduced to study wave turbulence 0], but 
these studies only lead to weak turbulence. Discrete lattice are attractive for 
two reasons. First they are simpler since they contain only a limited number 
of degrees of freedom (and therefore modes), but, more importantly because 
they have exact localized solutions which could play the role of vortices, the 
discrete breathers opening the possibility to study an analogous of fully 
developed turbulence. 

Many interesting experimental results on turbulence have been obtained 
with a simple model flow, the Von Karman flow produced in the gap between 
coaxial disks [^. We investigate here a nonlinear one-dimensional lattice 
excited by both ends, in a configuration which drives a large scale structure 
in the system, similar to the vortex of the Von Karman flow and we show 
that it exhibits remarkable similarities with turbulence, both for local and 
global properties. 

2 The one-dimensional lattice model. 

The system that we study is a one-dimensional chain of anharmonic oscilla- 
tors, coupled to each other by harmonic interactions. We have chosen Morse 
oscillators because, the presence of a cubic term in its expansion leads to a 
suitable nonlinear coupling between the modes, and moreover the lattice of 
Morse oscillators has localized oscillatory solutions which have interesting 
properties In analogy with the Von Karman flow, we want to drive a 
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finite lattice from botli ends, i.e. to inject energy in tlie system tlirougli an 
external driving. As in actual turbulence experiment where the viscosity of 
the fluids contributes do dissipate energy, we also need to introduce a dissi- 
pation mechanism so that the system can reach a steady state. In analogy 
with the dissipation term vV'^v of the Navier Stokes equations, where v is 
the kinematic viscosity, we introduce a dissipation term proportional to the 
second order finite difference of the velocity of the oscillators, so that the 
equations of motion of the nonlinear lattice are 

jjn - K{yn+1 + Vn-l - 2yn) - 2^06"^" (e"^" - l) - Z^(yn+l + Vn-l - ^i/n) = , 

(1) 

where y„ denotes the position of the n^^ oscillator, K is the coupling constant 
between neighboring sites, loq determines the depth of the Morse on-site 
potential, and u isa coefficient which plays the role of the kinematic viscosity. 
As usual, the dot designates a time derivative. 

Similarly to the Von Karman flow where the medium is excited from 
both ends by a set of rotating disks which create a large scale structure in 
the system, we consider a finite lattice of + 2 particles and we excite it by 
imposing the displacements of the two end particles, particle and particle 
A^ + 1. Their motion is chosen according to a linear mode of the lattice, with 
an amplitude A, a wavevector k = p7r/{N + 3), where l<p<A^-|-lisan 
integer, and a frequency co = \J 2u;q -|- sin k/2 (assuming that particles 
— 1 and A^ + 2 are fixed). Therefore the driving of the system is defined 
by two parameters, "p which fixes its wavevector and frequency and A which 
determines its amplitude. All simulations start from a lattice at rest and the 
equations of motion are integrated with a fourth-order Runge-Kutta scheme. 

Indeed, for such a nonlinear lattice, the fact that the driving corresponds 
to a single mode of the system does not imply that its response will be 
merely a single mode. Similarly to the Von Karman flow, driven by two 
co-rotating disks, which has a very complex response, the energy density 
nonlinear lattice also shows a complex spatio-temporal pattern, as displayed 
in Fig. Localized modes that move in the system, interact with each 
other, and die out due to dissipation are observed, and they appear to 
play a role very similar to the vortex filaments in fluids. The point of this 
study is to observe and analyze the response of the nonlinear lattice to 
a simple excitation. More complex cases could easily be imagined, such 
as an excitation with two different frequencies at both ends, which would 
be similar to some fluid experiments where the two disks have different 
velocities, but, as shown below the one-frequency driving is sufficient to 
yield interesting and puzzling results. 
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Figure 1: Energy density in a = 256-particle lattice excited by a single 
mode with p = 8, A = 0.12. The parameters of the model are loq = 1, 
K = 1, v = 5 10~^. The spatial dependence of the energy per unit celle 
e„(t) is shown with a grey scale (white corresponding to and black to a 
high energy density). The horizontal axis corresponds to time and extends 
over a time interval of 20 480 t.u., and the vertical axis gives the position n 
in the lattice. 



3 Local analysis. 

Numerical investigations of spatial the Fourier spectrum of the energy den- 
sity e{K,t) in this simple model have exhibited properties analogous to tur- 
bulence, such as a power law spectrum e(K) = Eqk°^ with a ~ 2.2 which is 
remarkably insensitive to the parameters of the model or of the excitation 



in a very broad range of parameters [11|. When the excitation amplitude in 



creases even more, a transition to a spectrum decaying exponentially versus 
K is found when the large scale structure created by the excitation breaks 
down into many smaller scale localized modes, as shown in Fig. |l[ This is 
the regime that we study in this letter. In actual turbulence experiments the 
energy density cannot be directly measured, but the local velocity is easily 
accessible. A map of the velocity field would require an array of probes that 
perturb the flow, and this is why the investigations only evaluate velocity dif- 
ferences measured in 2 points separated by a distance r, which provide some 
information on the spatial structure of the flow, and are also related to the 
viscous dissipation. For the nonlinear lattice it is easy to calculate an equiv- 
alent quantity in the numerical simulations u{5n,t) = yno+Snit) — yno(*)) 
where no is a fixed reference site and 5n is an index that measures the dis- 
tance between the two points, in units of the lattice spacing. It should be 
stressed that this quantity is local because there is no spatial averaging over 
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Figure 2: (a) Probability density functions of the velocity differences for 
6n = 2 (top plot) to 5n = 40 (bottom plot), by steps of 6n = 4. The points 
correspond to the histograms obtained by numerical simulations and the full 
lines are the fits of P{u) according to Eq. (|2|). The parameters of the model 
are = 256, loq = 1, K = 1, v = 10~^ and A = 0.12. The reference site 
is no = 57. For better lisibility, each distribution is shifter by -0.5 units, 
(b) Variations of the fitting parameter q versus dn/v for different viscosity 
parameters {u = 5 10~^, circles u = 1 10~^, stars, v = 5 10~^, x, = 2 10~^, 
squares) 



the reference site position uq. 

In a typical numerical experiment, we simulate the dynamics of the sys- 
tem for a time interval of 2 10^ time units (t.u.). The first 5 10^ t.u. are 
discarded in order to allow the system to reach a steady state before starting 
the recording of u{5n, t). Then we store u{dn, ti) at times ti separated by an 
interval tj+i — tj = 40 t.u. which is large enough to ensure that successive 
values are statistically independent. This gives us a set of 375000 values Ui 
for each 5n. We evaluate the reduced probability density of this distribution 
by first calculating ttj = (uj — {ui))/a where (•) designates a time average 

and o" = (uf) — {ui)"^ is the variance of the distribution. Then we make 
a normalized histogram of the ■Uj values in 100 bins to get the normalized 
probability density P{u). 

Figure ^a shows such histograms for 5n varying from 2 to 40. The 
analogy with similar figures deduced form experiments in fluid turbulence 
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is striking (see for instance Fig. 1 in Ref. |j^]). The same scenario from 
almost Gaussian distribution at large distances to stretched distributions at 
small distances is found. The analogy can be made more quantitative by 
analyzing the shape of the histograms. The curves of Fig. |a are very well 
fitted by the probability distribution deduced from nonextensive statistical 
mechanics [§, 12 1 



with P = 2/(5 — 3q), determined by the condition that the reduced probabil- 
ity distribution should have variance 1. As is a normalization constant, 
it should be noticed that each curve is fitted with only one free parameter, 
the parameter q > 1, which, in the formalism of nonextensive statistical 
mechanics measures the degree of nonextensivity. With g — > 1 ordinary sta- 
tistical mechanics is recovered and the probability distribution tends to a 
Gaussian. In the case of turbulence, a formula that generalizes @ has been 
derived to take into account the slight skewness of the distributions ||2[ . The 
distributions obtained for the nonlinear lattice also have a slight skewness 
so that the extended formula improves the fit. However, as it is already very 
good with Eq. @ we sticked to the expression that minimizes the number 
of parameters. 

The variation of q versus 6n is shown in Fig. ^3 for two values of the 
viscosity v and a mean amplitude (y) ~ 1.1 10~^ in both cases. As for 
the shape of the probability densities, the similarity with the experimental 
results on turbulence is strong. At small distances q is significantly larger 
than 1 (q 1.3 in our case) and it decays toward a value close to 1, but 
above 1, at large distances {q ~ 1.08 in our case). In turbulence experiments 
the variation of q versus distance is essentially determined by the rescaled 
distance r/rj, where rj is the Kolmogorov scale which is the typical size of the 
smallest vortices. In the discrete nonlinear lattice, the size of the smallest 
breathers is of the order of the lattice spacing, and, in the low viscosity 
range (high "Reynolds number") that we are investigating here, it does not 
vary significantly with the viscosity. This is why q only depends on 5n for 
the lattice as shown in Fig. B). 



4 Global properties of the dynamics. 

Our aim in this section is to examine to what extend the one-dimensional 
nonlinear lattice that we have investigated also shows the same "univer- 
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Figure 3: a) Probability distribution function (PDF) of the dissi- 
pated power Pdiss for different system sizes: from top to bottom N = 
256,512, 1024,2048,4096 The parameters of the model are lo^ = 1, K = 1, 
v = 5 10~^. The points correspond to histograms deduced form the nu- 
merical simulations and the full lines are fits with the generalized Gumbell 
distribution. For a better lisibility, each histogram is shifted by -0.5 unit 
along the vertical axis, b) Variation of the parameter b of the generalized 
Gumbell distribution versus the system size A'^ in a log-log scale. The error 
bars correspond to extreme values obtained in different simulations. 

sal" statistical properties that where observed in critical magnetic systems, 
self-organized criticality or turbulent flows when one investigate global quan- 
tities. 

In actual turbulence experiments, a global quantity which can be easily 
evaluated is the power injected by the engines driving the flow, for instance 
through the disks in the Von Karman experiment. In the numerical sim- 
ulations we have however the freedom to evaluate other global quantities 
directly in the bulk of the system. The power dissipated in the bulk 

n=3N/4 

Pdiss = XI ^'"^^'> ^ + 2/n-l - 2yn) Vn (3) 

n=N/A+l 

(where Fy is the viscous force) is particularly interesting because it probes 
the dissipation, an essential physical process for turbulence, at a global scale 
without being affected by the boundary effects. To avoid any perturbation 
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due to boundaries, we evaluate Pdiss by summing all the dissipative contri- 
bution from n = N/A + 1 to SN/i, i.e. in the central half of the system. 

As Pdiss is a global quantity, one could think that it would be smooth 
and featureless due to averaging. This is not the case. In addition to small 
amplitude fluctuations, the signal shows sharp peaks corresponding to strong 
bursts of dissipation. This feature shows up on the normalized probability 
distribution function (PDF) obtained from an histogram of the recorded 
values of -Pdiss in numerical simulations. Figure ^ show such histograms 
plotted in logarithmic scales as a function of = (Pdiss — (Pdiss)) /(^ where 
(Pdiss) is the time average of the dissipated power and Cdiss its standard 
deviation. The dissipation bursts give rise to a highly asymmetric PDF 
with a large tail on the negative side. The dissipation bursts are associated 
to collisions of the localized excitations, the discrete breathers, that are 
visible on Fig. |l[ Using an approximate expression of the breathers we can 
evaluate, in the continuum limit, the time average of the power dissipated 
by a breather. It shows that the dissipation grows with the third power of 
the breather amplitude, so that the large peaks observed during breather 
collisions generate strong bursts of dissipation. 

Figure |^ showing the space-time evolution of the energy density illus- 
trates a typical situation where one observes localized structures having a 
long life-time and a complex pattern of interactions. They travel on long 
distances introducing long range correlations in the dynamics of the lattice. 
In the context of the general studies of correlated systems which have been 
performed recently 0, systems showing such correlations have been called 
"inertial systems" and they display a strongly non Gaussian probability dis- 
tribution function of some global quantities. As shown by Fig. ^ this is also 
true for the nonlinear lattice that we investigate here. The strong asym- 
metry of the PDFs is similar to the observations made on turbulence and 
critical phenomena fs], ^. In order to evaluate this quantitatively, we have 
fitted the numerical PDFs with a generalized Gumbell distribution 

„b0 



U{e) =we^p ^a[be-e'"'j^ , (4) 

where w is the normalization factor and a, b are parameters of the fit. 
Figure H] shows that this expression gives a very good fit of the observed 
distributions, and is able to describe properly their evolution as a function 
of the system size N. In the limit 6 — > the distribution (^ tends to a 
Gaussian distribution ll{0) = exp[— 6*^/(2(70)] with a standard deviation 
(To = l/(a6^). The value 6=1 corresponds to the standard Gumbell distri- 
bution observed in a turbulent flow 0|, and larger values of b correspond to 



8 



even larger deviations from Gaussianity. Figure ^ shows that, even for large 
systems (A^ = 4096) the fluctuations of Pdiss stay highly non Gaussian. 

The origin of this non-Gaussian behavior can be understood from the 
same argument as for the case of turbulence or critical phenomena Q : there 
are long range correlations in the system which cannot be divided into sta- 
tistically independent regions. However the origin of these long range cor- 
relations appears to be different. While in the vicinity of a phase transition 
correlations are due to the existence of structures at all scales, in the present 
case, we observe small scale structures (localized modes) which can move on 
long distances and preserve their identities in collisions (see Fig. 0). How- 
ever the localized modes are not solitons and are therefore perturbed by the 
collisions, which limits the correlations to a finite range. As a result the 
PDF of the dissipated power depends on the size of the system and on the 
density of localized modes, both of which determine the average number of 
collisions that a localized excitation experiences when it travels inside the 
part of the lattice in which we compute the dissipated power (i.e. along N/2 
sites). 

The effect of the size of the system is shown in Fig. ^ which shows 
that b decreases with approximately with a power law b oc N~^'^''^ for 
128 < < 4096 so that the approach of a Gaussian distribution for large 
systems is very slow. The study of the effect of the mean amplitude of the 
displacement, which is correlated to the density of localized modes in the 
system, shows that a strong non-Gaussianity of the power dissipated in the 
system persists for variations of (y) that extend on more than two orders of 
magnitude. 

Consequently, contrary to the case of turbulence where a universal form 
of the fluctuations of a global quantity, independent of the Reynolds number, 
was found Q, the shape of the PDF of the power dissipated in the nonlinear 
lattice changes with the size of the system or amplitude of the excitation. 
However it stays strongly non- Gaussian and well described by a generalized 
Gumbell distribution for all the cases that we could investigate. One may 
wonder what are the differences between the turbulent flow and our model 
system. One possibility is the existence in the turbulent flow of structures at 
all scales, including the integral scale, i.e. flow patterns that have a charac- 
teristic size equal to the size of the system, which is not true for the nonlinear 
lattice. One may ask however whether larger Reynolds number would not 
break the large vortex into smaller filaments, leading to a gradual change 
in the shape of the PDF. Such a regime, if it exists, is probably beyond 
experimental possibilities, while, for the one-dimensional nonlinear lattice, 
the simulation can examine a very broad range of dynamical behaviors. 
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5 Conclusion 



In this work we examined a simple model related to turbulence. The statis- 
tics of its local and global properties are remarkably similar to the experi- 
mental results of actual turbulence so that it could be useful for its under- 
standing because it is sufficiently simple to allow detailed studies and test 
ideas. The remarkable features of the local and global properties could per- 
haps be understood in a single framework. Global properties of turbulence 
have been analyzed by Bramwell et al. 0] by pointing out that the flow can- 
not be divided into mesoscopic regions which are statistically independent. 
Such a "long range interaction" is precisely what leads to non-extensive 
thermodynamics |13] which models so well the local properties. 
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